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Introduction 

The homogeneous electron gas is a model which allows to study pure electron correlation 
without any interference with the multiple-scattering problem of real molecules, clusters, 
and solids. |^ The quantum kinematics of this phenomenon is hidden in the reduced 
density matrices and their diagonals. ||, ||, ^ The simplest reduced densites are (only 
the spin-unpolarized ground state is considered here) the 1-body momentum distribution 
n{k), recently parametrized in terms of the convex Kulik function, ^ and the 2-body 
quantities (711 (r) and being the non-negative pair densities (PDs) for electron pairs 

with parallel respectively antiparallel spins and with an interelectronic distance r. [0, ITUI, ITT 



All these quantities depend parametrically on the electron density 3/(47rrg). Correlation 
induced properties of n{k) (normalized as {'^/N)Y2y.n{k) = 1) are: (i) its nonidempotency 
< n{k) < 1, measured by a quantity c = {2/N) '^]^n{k)[l — n{k)] < 1 (which is christened 
here Lowdin parameter |jl2|) and (ii) the quasi-particle weight zp = "^,(1") — 77,(1"'") < 1 
(fc is measured in units of the Fermi wave length l/(ars),a = (4/97r)^/'^). Correlation 
induced properties of the PDs with their asymptotics (7||,|j(oo) = 1 are besides the oscillatory 
behavior for r 00, for small r the Fermi hole g^iir) < 1 with its on-top properties 
(yf|l(0) = 0, 5''|(0) = 0, and a characteristic curvature g'^}^{0) > and the Coulomb hole 
g-\l{r) < 1 with its characteristic on-top value g^-\{0) < 1. The on-top curvature of the Fermi 
hole is a local measure of the correlation strength. |T3|, [1^, |T6|, |T^ With the spin-summed 



PD g = [(yfji + 5'tj,]/2, particle-number fluctuations in spatial parts (fragments, domains) of 
the system can be discussed as another correlation index with the conclusion "correlation 



suppresses fluctuations" . [T9|, ^ |2T|, ^ |2J] For the ideal Fermi gas (r^ = 0) it is 



n(0)(A;) = 1 - 0(A;), = 0, = 1, and (/{^(r) = 1. 

The virial theorem provides a relation between the kinetic and the interaction energy, 
which follow from n{k) and g{r), respectively, giving thus for their r<. dependence an integral 



relation. p6|] In this paper another relation between n(k) and g{r) is derived. This deriva- 
tion is based on the successful parametrization of the PDs in terms of Overhauser 2-body 
wave functions (geminals), which are the scattering state solutions of an effective 2-body 



Schrodinger equation with an appropriately screened Coulomb repulsion, [g^, t3> HI' S 03 
It is furthermore based on the assumption that these Overhauser geminals can be used to 
represent also the 2-body reduced density matrix (2-matrix) 72(1|1', 2|2'), the digonal of 
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which gives the PD. Now the idea was to obtain the 1-matrix 71 from the contraction of the 
2-matrix 72. Unfortunately for extended systems this does not work, because this contrac- 
tion is not size-extensive. This problem is easy to overcome with the help of the cumulant 
, ^ ^ It defines by 72 = A7171 — x (A = antisymmetrizer) the 



expansion. |20|, 

size-extensively normalizable and contractable cumulant 2-matrix x = ^7i7i ~ I2, which is 
here thus represented in terms of Bessel functions (from ^47171) and Overhauser geminals 
(from 72). From the normalization of x follow sum rules (SRs) for the scattering phase shifts 



of the Overhauser geminals. |31| These (Friedel like) SRs are generalized in this paper by 
using the contraction properties of x, such that (at least in principle) n{k) can be calculated 
supposed the phase shifts of the Overhauser geminals are known. 



Basic notation and normalization sum rules 



The singlet(+ for even /)/triplet (— for odd /) components of the PD in terms of Over- 
hauser geminals Ri{r,k) are 



The fc-average is defined by 



<■■■>= 



N 



2 Qd^k 

Njo Wr 



{I, mi). 



d{k^ 



(1) 



(2) 



with the normalization volume Q and the density g = N/Q = I/Stt^. Wave lengths are 
measured in units of the Fermi wave length kp = l/{ars), a = (4/(97r))^/^. In this paper 
the Overhauser occupancy is slightly generalized by 



m = -Tf^iK,k), MK,k) = n( 



-K + k 

2 



-K-k 

2 



(3) 



N 

because it arises here not from the idempotent n^^\k) = 0(1 — k) of the ideal Fermi gas, 
but from the nonidempotent n{k) of the interacting electron gas. 0, H It is /i(0) = 8(1 — c). 
The -RL(r, k) = Ri{r, k)YL{er) are the (scattering state) solutions of the 2-body Schrodinger 
equation (the center-of-mass motion separates completely) 



[-A + v±{r)-e]RL{r,k)=0, v±{oo) = 



(4) 



with an appropriate repulsive interaction potential v±{r) = — + ■ ■ ■ , possibly different for 



" + " or " — " . |2^, 0, ^ The success of the Overhauser approach means, that 
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local interaction potentials f+(r) and f-(r), if not being exact, are at least reasonable 
approximat ions . 

One generalization of Eq. (|l]) concerns inhomogeneous systems. Here another gener- 
alization is considered, namely the representation of the 2-matrix in terms of Overhauser 
geminals: 

7±(R|R',r|r') = 0'WJ2% < ~^LL'{R~n',k)RL{r,k)Rl,{r',k) > (5) 

R = |(ri + r2) is the center-of-mass coordinate and r = ri — r2 is the relative coordinate. As 
/i(/c), the /-independent but /c-dependent weight in Eq.(|l]), also here the occupancy matrix 

follows with Eq. (j^) from the momentum distribution n{k). 

Arguments in favour of Eq.(|^) are: (i) The diagonal elements give the Overhauser PD 
Eq.(|I]), 7(R|R, r|r) = Q'^g±{r), and (ii) in the cumulant partitioning of 7± the generalized 
HP term is given by the same expression as Eq.(^ with only RL{r,k) replaced by 
jhikr) = ji{kr)YL{er) because the natural orbital are plane waves and using the L-expansion 
of a plane wave, exp(ikr) = 47r i'ji(A;r)F£(efc). 

So the cumulant 2-matrix x± = ~(7± — 7±^) is given by 



X±(R|R',r|r') = -^,2(4^)2^ < fiLAK-R\k)[R^{r,k)RUr\k) - j^{kr)jUkr')] > 



^L,L' 

(7) 

Por the normalization 



J rfVx±(R|R,r|r) = -g^ J d'r J^^l^ < fi{k)[R^{r, k) - jf{kr)] > 

-J2r / dk fi{k)7][{k) = ±gc, (8) 



= Q 

cf. Ref. 1^ 



Contraction sum rules 



With the contracted cumulant 2-matrices and their Pourier transform 

X±(|ri -r;|) = j d\, x±(R|R',r|r')[,^^ , = ^fP^ d\ e-"^^x±ir), (9) 



r2<-R— *oo 
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the contraction SRs are 

X±{k) = T Qn{K)[l-n{K)]. (10) 

They follow from 7±(k) = g n{K)[^N ± 1], 7^^(ft) = g n{K)[^N ±n{K)], and x± = — 7±, 
showing explicitly the cancellation of the non-size extensive terms. Eq. ([10|) says that n{k) 
can be calculated, if the Ihs is known. With x±(0) = jj J2kX±{'^) = TQC the special phase 
shift SRs of Ref. are contained in the more general SRs (|T0|) as special cases. 
With the identity (following from Eq. (H) and generalizing Eq. (13) of Ref. ||31||) 



d_ d_ 
dr dr' 



dr 



-[f±(r) - f±(r')] 



dr' 



(11) 
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(where i?^, = Ry{r',k)^ Rl= dRL{r,k)/dk'^, and h.c. means complex conjugate together 
with an exchange in both coordinates and indices) from Eqs.(P) and (P) it turns out 

X±ir) = xi{r)+xl{r) = -^'^^Y.^ < (^^' ^) i^LL'ir,k) + i?LL'(r, A;)] > . (12) 
The A and B matrices are defined by 

ALL'{r,k) + BLL'ir,k)=An J rfV [R^{r + r' , k)Rl,{r' , k) - j^{kir + r'))jl,{kr')] (13) 

where An' results from the 1st term of Eq. ([TI|) and -Bll' from the 2nd one. With the Gauss 
theorem it is 



ALL'{r,k) = 6ll'^^R^ 



dRi • <9 • Oil • d 

-g^Ri -R^d^R^-dRJi 



(14) 



(where Ri = Ri{R,k) and the r-dependence disappears) and with the analysis of Ref. ||3T 
it is finally 

(15) 



ALL'ir, k) = 6LuAi{k), Ai{k) = '^v'lik). 



So 



X±W = VEt|E^''''^'<M^,^)A(^)>, KK,k)= /"^/i(K,k). (16) 

K J 

Because of x±(0) = and Eq. (|^), it results x±(0) = —g'^ XIl < l^{k)Ai(k) >, in agreement 
with the (Friedel like) phase shift SRs of Ref. ||3l|, saying x±(0) = TQC. The Fourier 
transform 

xii^^) = -g' < 2'-"(2«^' >, (17) 
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enters (together with the B term) the more general SR (plQl). 

As aheady mentioned, the B matrix of Eq. (|I3]) follows from the 2nd term of Eq. (|Tl|) 



BLL'ir, k) = -An j d\' ^ K(|r + r'|) - v±{r')] x 
R^ir + r', k) Rl, (r', k)- R^ (r + r', k)RUv', k) 



So 



and 



Xl{r) = -q" 



(18) 



(19) 



(20) 



2" """ ^L,L' N ^1^.1 

Note X]kX±('*) = 0- Thus the B term does not contribute to the normalization of n{k) as 
it does not contribute to the normalization of the PDs. 

The Fourier transformed interaction potential and Overhauser geminals 



v^{k) = q d^r e-^'^''t;±(r), Rl{i^, k) = q d^r e-''"'"i?L(r, k) 



(21) 



determine the Fourier transformed B matrix 
in 2 



BiL'ii^i^k) = - 



Rl{k,2, k) Rl, (ki, k)- R^ (ki, k)Rl,{K,2, k) 



(22) 



1^12 = l*^! — ^^^21- Eq.([2^) has to be inserted into Eq. (pO|) . Then with (^12 = 6162) 

Pl"{C,\2)v±{l^\2) 

and with 

\x{2{k - Ki), k) = Pi^ ieiek)fii^ (k, ki, k) 
(the invariance of /i(K, k) by the replacement k — k makes h even) and with 



Cii^i' 



dn m)p,,{OP,{C) 



(23) 



(24) 



(25) 



(because / and /' have the same parity, these coefficients vanish for odd h) an expression 
arises, 

v^{k,i, K2)Ri{k2, k) Rl' (ki, k)- Ri [ni, k)Ri>{K,2, k)v^{Ki, K2) >, (26) 
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which vanishes, because the L-L' double sum runs over an antisymmetric matrix (each 
element 1,1' is compensated by its mirror party). Thus, it finally results 

< ^^^%^4{k) >= ±n{k)[l-n{k)]. (27) 

The occupancy weight in Eq. ( P?] ) can be written in terms of n{k) as 

fi{2t,,k) = j ^n(|Av + k|)n(|Av-k|), (28) 

or with n{k) = |p/ d^r '^'"^^ /(r) in terms of the dimensionless 1-matrix /(r) as 

1 9 /" ,•:! ^■x sin kIti + rol sin fclri — r^l x x . ^ 

k) = -/ j Sndh, k\^^_,p ^^^)f^^-)^ (29) 

For the momentum distribution the parametrization in terms of the Kulik function can be 
used. [0, i 



Conclusions 



The contraction SRs (^) are nonlinear integral equations for n{k) supposed the phase 
shifts of the Overhauser geminals Ri{r, k) are known from the effective interaction potential 
v±{r) = ^ + ■ ■ ■ . These integral equations have to be solved selfconsistently. Because the 
PD is (via the Ri{r, k)) a functional of f±(r) and n{k) is a functional of the Ri{r, k), thus 
the total energy becomes a functional of v±{r). 



The success of the Overhauser approach seems to confirm the possibility of a pair-density 



functional theory, the idea of which is presented in Ref. |3g]. Whereas therein the geminal 
occupancy was erroneously assumed to be 1 and according to occupied and unoccupied 
geminals respectively (analog to the aufbau principle of the Hartree-Fock and the Kohn- 
Sham schemes), in the Overhauser approach a geminal occupancy fi{k) is used being quite 
different from a step function. For PD or 2-matrix geminals there is no aufbau principle. 
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